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ABSTRACT 

Recently, low light level charge coupled devices (L3CCDs) capable of on-chip gain 
have been developed, leading to sub-electron effective readout noise, allowing for the 
detection of single photon events. Optical interferometry usually requires the detec- 
tion of faint signals at high speed and so L3CCDs are an obvious choice for these 
applications. Here we analyse the effect that using an L3CCD has on visibility pa- 
rameter estimation (amplitude and triple product phase), including situations where 
the L3CCD raw output is processed in an attempt to reduce the effect of stochastic 
multiplication noise introduced by the on-chip gain process. We establish that under 
most conditions, fringe parameters are estimated accurately, whilst at low light levels, 
a bias correction which we determine here, may need to be applied to the estimate of 
fringe visibility amplitude. These results show that L3CCDs are potentially excellent 
detectors for astronomical interferometry at optical wavelengths. 

Key words: instrumentation: detectors - instrumentation: interferometers - tech- 
niques: interferometric - methods: statistical - methods: numerical. 



1 INTRODUCTION 

Optical interferometry requires the detection of an interfer- 
ence fringe pattern generated when starlight collected by 
separated apertures is recombined with a spatially or tem- 
porally modulated optical path difference (OPD). If the to- 
tal distance travelled by the two beams of light is equal, 
or the OPD is a whole number of wavelengths, constructive 
interference will result in an intensity maxima. Likewise, if 
the OPD is an odd number of half wavelengths, destructive 
interference results in an intensity minima. In a typical in- 
terferometer, intensities are measured for a range of OPDs, 
with the OPD differing by a distance of several wavelengths 
about the zero OPD. Spatial and temporal modulation of 
the OPD can both be used by optical interferometers, for 
exampl e the Amber instrument on the VLTI (spatial mod- 
ulation, iR/jbbe^Dufx^ the COAST (tem- 
poral modulation. lBaloVmeta!T Tl999)h An oscillating sig- 
nal or fringe pattern results, the contrast and phase of these 
fringes providing information about the source, including the 
source diameter and shape, and, when the light is combined 
using three or more apertures, information about asymme- 
try within the object. The entire fringe must be sampled 
within a few times the time-scale of atmospheric perturba- 
tions (typically 10 ms in the optical regime), resulting in 
low signal intensities, and for all but the brightest stars the 
signal is very faint. To date, the faintest reported optical 
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measurement using a separate element interferometer is an 
objec t with visible magnitude of about 8.5 (I magnitude 6.1) 
iHaniff et al.ll200(J) . Increasing the sensitivity of interferom- 
eters is therefore essential if they are to be used to observe a 
wide range of scientific sources. One way of doing this is to 
observe over a wide range of wavelengths, allowing a greater 
light throughput. However, to maintain a useful coherence 
length, the bandwidth of light on any detector element must 
be kept small, requiring independent detection of many dif- 
ferent wavelengths (spectroscopic detection). This will re- 
quire at least a one dimensional array of detector elements, 
which may be expensive and fragile if conventional photon 
counting detectors (APDs and PMTs) are used. 

CCDs ought to be ideal detectors for spectroscopic in- 
terferometry because they are stable devices available in 
large-format arrays, allowing one CCD to replace many sin- 
gle element detectors. However, their major shortcoming is 
readout noise, i.e. the additional signal added at the on-chip 
output amplifier where the photo-electrons are counted. A 
typical interferometer can require CCD readout rates of up 
to and above 1 MHz since there will be many spectral chan- 
nels, and each of these must be sampled many times within 
an atmospheric coherence time. Currently, the best noise 
level achieved usin g standard CCDs at these readout rates 
is typically 10 e~ llJerram et aljl200lf) . Signal levels lower 
than this are then swamped by noise. 

Low light level CCDs (L3CCDs) provide a solution to 
this problem, with an on-chip gain ampl ifying the signa l 
prior to the on-chip readout amplifier llJerram et alJl200lfl . 
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resulting in a sub-electron effective readout noise. This, com- 
bined with the high quantum efficiency (QE) of L3CCDs (up 
to 90 percent), which comp ares favourably with the QE of 
APDs IITakeuchi et al.lll999T) . makes them prime candidates 
as interferometric detectors. Being array detectors they can 
be used for spectroscopic detection allowing greater use of 
the available light, but they do have the disadvantage of in- 
troducing an additional source of noise due to the stochastic 
multipl ication process. How ever for one application, we have 
shown feasden et al]l2003l) that the magnitude of this noise 
can be reduced if the L3CCD output signal is treated in the 
correct way, by thresholding. 

In this paper we investigate the use of L3CCDs for inter- 
ferometric fringe detection. In most situations we find that 
no bias correction is required to estimate fringe parameters 
correctly. In some cases, a correction is necessary when es- 
timating visibility amplitude but we provide this correction 
here. We also develop a suitable treatment allowing unbiased 
estimation of bispectrum (closure) phase. 

The layout of the paper is as follows. In section 2 we 
present a background to interferometric fringes and the use 
of L3CCDs, as well as our modelling techniques. In section 
3 we present results and our conclusions are summarized in 
section 4. 



2 INTERFERENCE AND L3CCDS 

2.1 Introduction to interferometric fringe 
production 

In interferometry, the raw image we obtain is an interfero- 
metric fringe pattern, which will show a periodic intensity 
variation due to either spatial or temporal OPD modula- 
tion depending on the type of interferometer. An expression 
for a g eneral interferometric fringe pattern is given bv lScottJ 
il997t) (Fig. although for the purposes of this paper we 
use a simplified version ignoring atmospheric effects, and 
assuming a monochromatic signal (Eq. . 

Once a fringe pattern has been recorded, it is useful to 
obtain the visibility amplitude and (for space based interfer- 
ometers) visibility phase. If three or more separate apertures 
have been used, it is possible to measure the bispectrum (or 
closure) phase. 



2.1.1 Typical light levels 

The typical light levels required for measurement by an in- 
terferometer are very low, which is why conventional CCDs 
with higher readout noise are not acceptable as interfero- 
metric detectors. For example, observation of a sun-like star 
with a lossless interferometer at a distance of 1 kpc, with 
two telescopes of size 1 m diameter, a single pixel detector 
with a sampling rate of 5 kHz, with a bandpass of 30 nm 
centered at 700 nm, and a detector with 50 percent quan- 
tum efficiency would produce an average detected signal of 
about 0.7 photons per detector pixel per readout from the 
interference fringe. It is therefore essential that an interfer- 
ometric detector must be able to detect individual photons 
if anything other than bright sources are to be observed. 
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Figure 1. Example of a simplified fringe pattern with a maximum 
visibility amplitude of 0.5. The enveloping shape of this fringe 
pattern is due to a finite signal bandwidth, AA/A = 0.2 in this 
case. 



2.2 L3CCD multiplication noise 

The on-chip multiplication process in an L3CCD is stochas- 
tic and introduces extra noise into the output signal of the 
L3CCD before readout, in addition to the Poisson noise due 
to the detected photons. This reduces the signal-to-noise ra- 
tio (SNR) of the raw output by up to y2 . By processing the 
raw measured output feasden et all2o'()3T) we can reduce the 
effect of this stochastic noise at low light levels, the reduc- 
tion depending on the processing strategy that we choose. 
Increasing the SNR of the L3CCD output signal then allows 
us to improve our estimate of interferometric visibilities. 

We use the processing strategies suggested by 
iBasden et alj |2003) which were developed for photomet- 
ric applications, namely the "Analogue", "PC" and "PP" 
thresholding strategies which are defined as follows: 

Analogue: The L3CCD output is divided by the mean 
gain. 

Photon counting, PC: If the L3CCD output signal is 
above some noise threshold, it is assumed to represent one 
photon. 

Poisson probability, PP: Thresholds are placed at the 
positions where the L3CCD output probability distribution 
for a given mean input of n (integer) Poisson photons crosses 
with the output distribution for a mean input of n+1 Poisson 
photons. 

Additionally we here explore a further thresholding 
strategy denoted as "Uniform" , defined as follows: 

Uniform: Places the L3CCD output signal into evenly 
spaced thresholds, the spacing being dependent on the mean 
light level (estimated from the fringe), chosen such that on 
average, at a given light level, we will estimate the mean 
photon input correctly. This allows us to make better use of 
the fringe data. Further details of this thresholding strategy 
are given in appendix 1X1 
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2.3 Numerical models 

We investigate the effect of using an L3CCD both with and 
without processing of the raw output for a range of visi- 
bility amplitudes and light levels, categorizing the perfor- 
mance of L3CCDs as interferometric detectors, particularly 
in the estimation of visibility amplitude and bispectrum (clo- 
sure) phase. Any additional bias introduced to these mea- 
surements by the multiplication process needs to be deter- 
mined. In this paper we investigate the effect of L3CCDs 
using Monte-Carlo simulation as follows: 

(i) Multiple realizations of different fake interferometric 
fringes were generated. 

(ii) Poissonization of the fringe signal was simulated. 

(iii) The L3CCD stochastic multiplication process on the 
Poissonized signal was simulated. 

(iv) The simulated L3CCD output was processed using 
techniques known to reduce the noise introduced by the mul- 
tiplication process (i.e. Uniform, PP, PC). 

(v) The effect of these processing strategies on fringe vis- 
ibility estimation (visibility amplitude and phase, and bis- 
pectrum phase) was investigated using standard methods, 
and checked to see whether the estimates were the same as 
when using a perfect Poissonian detector. 

The interference fringe pattern is modelled assuming a 
two element interferometer, using the following expression, 



1e+07 - 



Ij = n [1 + V cos (2njm/M + </>)] 



(1) 



where 7j is the intensity of the fringe signal in the j th pixel, 
for a fringe with a mean intensity \i and visibility ampli- 
tude V, with phase <j>. The fringe extends over an array of 
size M pixels (M = 256 in our case) and contains m whole 
wavelengths within this array, so that we are not affected 
by aliasing when Fourier transforming. In our calculations, 
we ensured that there were always more than eight sam- 
ples between adjacent peaks of the fringe, so that they were 
well sampled. Poisson photons with mean intensity Ij were 
then produced from the theoretical fringe for a given OPD, 
giving the fringe as it would appear to a detector before mul- 
tiplication and readout. By modelling the behaviour of the 
L3CCD multi plication regis t er usi ng a Monte-Carlo simula- 
tion following iBasden et al.l 1120031) . the effect of L3CCDs on 
visibility parameter estimation was investigated. 



2.4 Visibility amplitude estimation 

It is well established jWalkup fc Goodman! Il973f) that the 
mean power spectrum (square of the Fourier transform) of 
an interference fringe averaged over many realizations will 
yield a DC peak with height N 2 (where TV is the mean 
number of detected photons in a single fringe), a peak at a 
position corresponding to the fringe modulation frequency, 
and a white noise background due to Poisson noise and 
L3CCD multiplication noise, as shown in Fig. [5] In the 
case of a two element interferometer considering only pho- 
ton noise, an unbiased esti mate of the visibility amplitude 
JWalkup fc Goodmar]|l973Tl is obtained using 
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Figure 2. Sample powerspectra of interference fringes produced 
from L3CCD simulations. The lower (black) trace is produced by 
a two element interferometer with a mean light level of 1 photon 
per pixel, and a visibility amplitude of 0.75 from the addition of 
20 fringe powerspectra. The upper (grey) trace is produced by a 
three element interferometer with a mean light level of one photon 
per pixel, and visibility amplitudes of 0.75, 0.5 and 0.25 from the 
addition of 200 fringe powerspectra. The white noise background 
can be seen clearly, and shows that the multiplication process 
does not alter the slope of this noise. 



where E[D (2) ] is the power spectrum of the fringe measured 
at the modulation frequency, u, or DC, and b is the noise 
background level. L3CCD multiplication noise may intro- 
duce a bias to the estimator depending on the output thresh- 
olding strategy and light level. 

In some cases, detector im perfections c an introduce a 
slope to the noise background iPerrirJ 2003). However, we 
have established that these effects are not relevant when 
using an L3CCD and so we do not consider such effects 
here; 

IWalkup fc Goodman] (11973ft give estimates for the num- 
ber of photon detections that are required to measure the 
visibility amplitude of a fringe to a specified accuracy. To 
measure a visibility amplitude V to a SNR of 1Z in the pho- 
ton limited case, we require 



N 



2TZ 2 
V 2 



(3) 



I VI = 2* 



1 E[DM (u)] - b 
E[D(V(0)]-b 



(2) 



total signal detections. We therefore generate fringes and 
sum the power spectra until the total number of Poisson 
photon detections allows us to reach the required SNR (in 
our case 50). We then use the summed power-spectrum to 
estimate the visibility amplitude. We investigate visibility 
amplitudes between 0.02 — 1 at mean light levels between 
0.01 — 20 photons per pixel, chosen in the range where 
IBasden et alJ (I2003T) find L3CCD output thresholding has 
a positive effect on the SNR. 



2.5 Visibility phase estimation 

The visibility phase is meaningful only for space based inter- 
ferometers, since atmospheric perturbations render it mean- 
ingless unless the phase introduced by the atmosphere is 
known. It is estimated from the phase of the Fourier trans- 
form at the fringe modulation peak. When averaging multi- 
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pie realizations of the fringe pattern, we average the complex 
value of the peak, not just the phase. This ensures that we 
give the correct weighting to each individual value, depen- 
dent on the number of photons detected. 



2.6 Bispectrum phase estimation 

When light from at least three apertures is combined at 
different (non-redundant) OPD modulation frequencies, we 
can use the detected signal to retrieve some of the phase 
information by summing together the phases, cancelling the 
relative phase introduced by the atmosphere. To model this 
we use a simpl i fied m odel of the generalized fringe pattern 
given bv lScofrJ (1199711 . The fringe power spectrum will con- 
tain three signal peaks and a DC peak as shown in Fig. 
(top trace). The bispectrum phase is estimated from the 
phase of the average complex bispectrum product, with a 
bias correction that must be applied at low light levels. A 
suitable correction to apply to the measured bispectrum in 
the case w hen detection is g overned solely by photon noise 
is given bv lWirnitzerl lll985l) : 



< Q (3) (u,v) > 



E [D (3) (u,v)] +dN (4) 
-c 2 E[\D (2) (u)\ 2 + \D (2) (v)\ 2 
+ \D (2) {-u-v)\ 2 } 



where ci = 2, oi = 1, < Q l ' (u,v) > is the unbiased esti- 
mate for the bispectrum and E [Z?' 3 ^ (u, v)~\ is the expected 
triple product of the fringe (the product of the Fourier trans- 
form of the data evaluated at the three signal peaks). We 
now evaluate the phase of (u,v) instead of the photon 
noise biased value of (u,v), with u and v correspond- 
ing to the different fringe frequencies, D^(u) representing 
the fringe power spectrum, and N being the total photon 
count. Eq.^Jshows that unlike the power spectrum, the pho- 
ton bias introduces both constant and frequency dependent 
terms into the bispectrum. 

If the form of the noise is not Poissonian, as is the case 
when using an L3CCD, a different bias c orrection may be 
requir ed. We are able to use results from iPehlemann et alJ 
( 1992) to determine the bias correction for the case when we 
use the raw L3CCD output, obtaining c\ = 6 and C2 = 2. 
The details, which involves considering the L3CCD output 
probability distribution, are discussed in appendix IU1 

If we do not know the detector statistics, it may still be 
possible to obtain an unbiased estimation of the bispectrum 
phase. To do this, we used the data to compute ci and C2 
in Eq. 2] This technique relies on Eq. 0] being the correct 
model for the form of the bias correction, and assumes that 
the power spectrum can be written as 



E[D (2) (u)] = c PS ,o < J {2) (u) > +c PS ,i 



(5) 



where cps,o and cps,i are scalar coefficients, and < 
J (2) (u) > is the average power spectrum of the normalized 
high-light-level (unbiased) fringe. 

The coefficients ci and C2 were obtained by evaluating 
the bispectrum at frequencies u, v chosen such that other 
terms disappear. We give further details in appen dix IO and 
such a scheme was initially proposed by IPehlemann et ail 
(1992). The form of the power spectrum assumed by this 
scheme is not always correct, though since differences will 



be small, we attempted to use this scheme to estimate the 
bispectrum phase. In most circumstances we find that a cor- 
rection can be found. 



2.7 Theoretically predicted visibility amplitude 
bias 

If we use PP or Uniform thresholding strategies on the 
L3CCD output we will introduce a systematic bias into our 
visibility amplitude estimations at some light levels since we 
place different weightings on fringe minima and maxima. We 
can get an understanding of how to treat this by consider- 
ing a theoretical model of the L3CCD out put probability 
distribution as given bv lBasden et al] 1.20031) : 
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_1 exp(-:r/g) 



n!(n - l)\g n 



(6) 



where g is the mean gain, /1 is the mean light level and x is 
the number of electrons at the L3CCD output after multi- 
plication. We can therefore find the estimated light level I in 
terms of the true light level fi for any thresholding strategy 
by writing 



a — 1 x — f a _ 



P(x) 



(7) 



where the a th threshold boundary is represented by f a and 
P(x) is given above. 

The classical definition for visibility is 

t r -/max -^min , \ 

= 1 — TT~ ' { ' 

at a given light level fi = (7 ma x + -Zmin) /2 and true visibility 
V. The bias in the estimated visibility amplitude can then be 
determined, and so corrected, giving us an unbiased estimate 
for visibility amplitude. We expect our visibility estimate V c 
to be equal to: 



K(M, V) 



with 



E a ,^S(x,n)(W + + W.) 



(9) 



S(x,n) — 



fi n x" x exp(— x/g) 
n\(n-l)\g n 



W± = exp(=FM0(l± V) n 
/max = V] aS(x,n)W+ 

a,x,n 

where the summation is over n = 1 — > 00, x = / a -i — ► 
/ a and finally a = 1 — * 00. We have used this theoretical 
model alongside our Monte-Carlo simulations to verify the 
predictions we make. 

In a similar way, we can compute the visibility ampli- 
tude bias when using a single threshold. Since we interpret 
any signal as a single photon, at a given true light level {i 
we measure the mean light level to be 1 — exp(— fi). The 
true light level at the fringe maxima and minima is given by 
I± — ± V), and so inserting this into the expression for 
visibility gives 



exp(~fi){exp(fiV) - exp(-^V)] 
— exp(— ^)[exp(pV) + exp(— fj,V)] 



(10) 
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where the symbols are as before. 



3 ACCURACY OF FRINGE PARAMETER 
ESTIMATION 

We separate our analysis into three light level regimes, which 
allows us to apply different L3CCD output thresholding 
strategies as appropriate for maximum reduction of the noise 
introduced by the L3CCD multiplication process. Our light 
level regimes are chosen according to: 

Low light levels: Much less than one photon per pixel per 
readout (< 0.1 photons per pixel). 

Intermediate light levels: Between about 0.1 — 20 photons 
per pixel per readout. 

High light levels: More than 20 photons per pixel. 

Read noise and gain were found to have little effect on 
the visibility estimations, provided the mean g ain was s ignif- 
icantly greater than the read noise as found bv lBasden et alJ 
(2003). If this was not the case, then the L3CCD output sig- 
nal was dominated by noise, leading to inaccurate visibility 
estimation. 

We find that all of the L3CCD output thresholding 
strategies leave the mean power spectrum of the fringe with 
a flat background. The height of this background is deter- 
mined by the variance of the signal and in a Poisson case it 
is equal to the number of photons detected (which is also the 
variance). When using an L3CCD analogue processing strat- 
egy, we find that the background level is twice the height 
of that from a pure photon noise case as expected since 
the signal variance is double that of the Poisson case. Sim- 
ilarly, with a single threshold we find that the background 
level is lower than in the pure photon case by an amount 
in agreement with the reduction in signal variance. Multiple 
thresholding strategies give a background level above the 
pure photon case and below the analogue case, the height 
depending on light level, in agr eement with the variance of 
the signal feasden et al J 12003ft which is dependent on the 
light level. 

3.1 Visibility amplitude estimation at low light 
levels (< 0.1 photons per pixel) 

As described bv lBasden et alJ d2003ft . if the mean light level 
is low (much less than one photon per pixel per readout) 
then it is most probably that either one or zero photons 
will land on a pixel during an exposure. At the L3CCD out- 
put we will then have either a relatively large signal corre- 
sponding to one photon, or a small signal, just due to noise 
from the electronic readout amplifier. We can therefore use a 
PC thresholding strategy treating every signal above a noise 
level as representing one photon and operating the L3CCD 
in photon counting mode. We find that the visibility ampli- 
tude estimations are unbiased for all visibilities provided the 
light level is kept low, as shown in Fig. [31 (d) . 

A uniform thresholding strategy will lead to the same 
results, since at low light levels it is virtually identical to a 
PC thresholding strategy. Similarly using an analogue strat- 
egy (the raw output), though requiring about twice as many 
detections t o achieve the same SNR due to a larger excess 
noise factor feasden et al$ 2003). will also provide unbiased 



visibility estimation. A PP thresholding strategy also pro- 
vides an unbiased estimate for visibility. 

As the light level increases, using a single threshold (PC 
strategy) causes us to lose estimation accuracy, since there 
is an increased probability of two or more photons being de- 
tected in any given pixel. These signals will be interpreted 
as a single photon event, and so the fringe maxima will be 
underestimated more than the fringe minima. Visibility am- 
plitude will be underestimated, and at high light levels will 
tend towards zero. Other thresholding strategies are then 
more appropriate. 

3.2 Visibility amplitude estimation at 

intermediate light levels (0.1 — 20 photons per 
pixel) 

At intermediate light levels, we cannot use a single threshold 
processing strategy as coincidence losses are large. However, 
we can still reduce the excess noise introduced by the multi- 
plication process by thresholding the L3CCD output signal. 

3.2.1 PP and uniform thresholding 

We can threshold the L3CCD output signal into prede- 
termined thresholds which are chosen to minimize the ex- 
cess noise factor introduced by the stochastic multiplica- 
tion process, for exa mple using the PP thresholding strategy 
( Basd en et aljl2003ft . Alternatively, we can use the raw data 
from a fringe to estimate the mean light level. We can then 
use this along with our knowledge of the mean gain to im- 
prove our input prediction, estimating the most likely value 
of the photon input by choosing appropriate thresholds tak- 
ing into account the mean light level (appendix^. For any 
given pixel, the number of photon events may be estimated 
incorrectly, though on average we will estimate correctly. At 
low signal levels, we will estimate correctly most of the time, 
minimizing the dispersion introduced by the multiplication 
process. As t he light level increa ses, we will lose accuracy 
progressively feasden et al.l 12003ft until there is no advan- 
tage in thresholding over using the raw L3CCD output. 

Using the PP or Uniform thresholding strategies at light 
levels between about 0.1 — 10 photons per pixel will intro- 
duce a small bias to the estimated visibility amplitude. We 
find that this bias is greatest at light levels of about two pho- 
tons per pixel, and that here the estimated visibility is just 
over 80 percent of the true visibility. This is independent of 
the mean gain, provided the gain is well above the readout 
noise. An explanation for this bias is given in appendix iBl 
and we find that our Monte-Carlo simulations agree with 
the theoretically prediced bias (section |2. 71 . We find that a 
simple bias correction can be applied according to 



y cor r/o v. 7 

v: /2 (i - h{n)) + h(p) 

where h(fi) = 1 — /j 3 ^ exp(— /x/3)/5 for a mean light level 
/j, photons per pixel, with V e representing the biased visi- 
bility estimated from the thresholded data, and V CO r being 
the corrected visibility amplitude. It is important to note 
that this is only an approximation generated by fitting to 
the data, though it allows the visibility amplitude to be esti- 
mated to about two percent accuracy, i.e. the visibility will 
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Figure 3. Plots showing visibility amplitude estimates when using an L3CCD as a function of light level (mean photons per pixel), for 
visibilities of 1, 0.75, 0.5, 0.2, 0.1, 0.05 and 0.02, with points representing Monte-Carlo simulation results (SNR 50), and curves being 
theoretical, (a) Results when using an analogue thresholding strategy. In this case no bias is evident, (b) and (c) Showing results when 
using a uniform thresholding strategy without and with a correction required to estimate amplitude correctly, (d) Showing results using 
a PC thresholding strategy. Here we see that this cannot be used at light levels greater than about 0.1 — 0.5 photons per pixel, (e) and 
(f) Showing visibility amplitude estimates when using a PP thresholding strategy, before and after correction for visibility amplitude 
bias respectively. 



be estimated to within (1 ± 0.02)Vtruc- A more accurate bias 
correction can be estimated by considering the L3CCD out- 
put probability distribution, though this requires analytical 
treatment to compute it. 

Monte-Carlo simulation results before the application of 
the simple bias correction, along with the variation of Eq. [5] 
and Eq. llUl with [i and V are shown in Fig.|H](b), (d) and (e) 
for the uniform, PC and PP thresholding strategies. Fig. [3] 
(c) and (f) show the visibility amplitude estimates after the 
simple bias correction has been applied. We find that the 
visibility bias is very similar for both the uniform and PP 
thresholding strategies. 



levels in all cases where there is ample SNR. Visibility am- 
plitude reproduction is shown in Fig. |3](a). 

3.3 Visibility amplitude estimation at high light 
levels 

At high light levels, there is no real advantage in threshold- 
ing the L3CCD output as this does not improve the SNR, 
and so an analogue processing strategy should be used, 
treating the L3CCD output as that from a conventional 
CCD. Visibilities (amplitude and phase) will be estimated 
with no additional bias. 



3.2.2 Raw output 

Using the analogue processing strategy provides an unbiased 
estimate for visibility, since we are not weighting the data 
in any way. This would seem to be a big advantage, except 
that the SNR can be reduced by a factor of up to v2 when 
compared with other thresholding strategies. We therefore 
recommend the use of this processing strategy at these light 



3.4 Phase 

The visibility phase estimate from a fringe produced by a 
two element interferometer is unbiased in all of the previ- 
ously mentioned readout modes (section 12.21 at any light 
level. It is only fringe peak height, and not position of the 
fringe that is altered by thresholding the raw L3CCD out- 
put, and so the L3CCD has no effect on the visibility phase 
estimation. 
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The spread in a p hase measurem ent estimated from 
many fringes is given bv lBuscherl (|l988) as: 



yTl < COS £ > 

where e is the deviation of the phase estimate from the 
true phase in each fringe pattern, and n is the number of 
phase measurements (fringes). We find that a generated 
by our simulati ons agrees well with the t heoretical SNR 
a" 1 = NV 2 /2 JWalkup fc Goodman! Il97l . when the ex- 
cess noise factor introduced by the L3CCD multiplication 
process after thresholding is taken into account. An exam- 
ple from our simulations is shown in Fig. 2] 

3.5 Bispectrum phase 

Our simulations show that when using an L3CCD we are 
able to estimate the bispectrum phase accurately at mean 
light levels down to about five photons per pixel when using 
the standard photon noise bias correction (Eq. |1J . At lower 
light levels, the bispectrum phase estimate appears to be 
biased (Fig.|SJ. If an analogue or PP thresholding strategy is 
used with the standard photon bias correction, we find that 
the estimate of phase tends towards zero as light decreases, 
even if photometric corrections are applied when using the 
PP strategy. 

In the analogue case, we are able to estimate phase cor- 
rectly by using the bia s correction that we adapted from 
iPehlemann et alJ <1992ft (Eq.gjwith ci = 6 and c 2 = 2). At 
very low light levels, a large number of fringe patterns need 
to be recorded to reduce the SNR to acceptable levels. 

Using a unif orm thresholding strategy and the standard 
IWirnitzerl il985T) bias correction, we find that the phase is 
overestimated slightly at low light levels as show n in Fig. |S| 
(a). A PC thresholding strategy used with the IWirnitzerl 
(1985) bias correction is found to overestimate the phase 
slightly up to light levels of about 0.2 photons per pixel, 
above which the estimate tends quickly towards ±7r as co- 
incidence losses become large. 

The failure of Eq. 0] when estimating the phase is due 
to the L3CCD output probability distribution no longer be- 
ing Poissonian, which is required for the correction to work 
correctly. However, by calculating appropriate values for the 
coefficients in Eq.0]suitable for an analogue processing strat- 
egy, we are able to obtain an unbiased estimation of bispec- 
trum phase when using the this processing strategy down to 
very low light levels. 

A similar correction for other L3CCD output processing 
strategies has not been determined due to the complicated 
nature of the probability distributions. However, by using a 
different approach to bispectrum estimation, using the dat a 
itself to calculate the bias terms llPehlemann et al] il992l) . 
section |2~B1 and appendix [CJ, we find that we are able to es- 
timate bispectrum phase accurately to light levels down to 
about 0.1 photons per pixel when using the PP or uniform 
thresholding strategies. This corresponds to the light levels 
at which L3CCDs are most likely to be used. Below this, 
errors in the individual bias terms, and invalid assumptions 
about the form of the power spectrum reduce the phase es- 
timate below the true value, as shown in Fig. |S] (b) . 

We recommend the use of the analogue processing strat- 
egy and our adapted bias correction as this gives an unbiased 



estimate of bispectrum phase. Alternatively, if the SNR is 
low or at very low light levels, we recommend the use of the 
Uniform thresholding strategy since this is able to reduce 
the effect of stochastic multiplication noise. The Uniform 
thresholding strategy will however lead to a slight overesti- 
mation of the bispectrum phase by between about 10 — 20 
percent at the lowest light levels. 



3.6 Uncertainties in gain 

A knowledge of the mean gain is required when using ei- 
ther the PP or uniform thresholding strategies. If this is 
not known precisely, threshold boundaries may be placed 
wrongly. Fortunately, due to the nature of visibility, being 
composed of both a difference and a ratio (Eq. [HJl, the effect 
on the amplitude estimation will be much less than the error 
in gain. 

We find that when using a PP thresholding strategy, 
the uncertainty in estimated visibility amplitude is inversely 
proportional to the uncertainty in the gain estimate. If we 
overestimate the gain by 20 percent, then in the worst case 
(for low visibilities (< 0.1) and light levels around two pho- 
tons per pixel), we will underestimate the visibility ampli- 
tude by at most four percent for a PP thresholding strategy, 
and five percent for a uniform thresholding strategy. For ex- 
ample, if the true mean gain was 1000, but we had estimated 
it to be 1200, and are trying to measure a true visibility am- 
plitude of 0.05 at a mean light level of 2 photons per pixel, 
we would estimate the visibility amplitude wrongly by about 
five percent, i.e. we would estimate the visibility amplitude 
to be 0.0475. 

For visibilities above about 0.1, and mean light levels 
fainter or brighter than two photons per pixel, the error 
in visibility is reduced. For PC and analogue thresholding 
strategies, uncertainties in the mean gain have no effect on 
visibility estimation. 

Typical uncertainties in the gain are of order one per- 
cent jMackav et alJl200lT) . and so in general there will be 
negligible uncertainty in the visibility estimate due to this 
uncertainty. Visibility phase estimates are unaffected by un- 
certainties in the mean gain. 



4 CONCLUSIONS 

We have investigated the application of L3CCDs to inter- 
ferometric fringe detection, and the way in which this af- 
fects visibility parameter esti mation. In doing thi s, we have 
introduced an additional (to [Basden et all (l2003f) l L3CCD 
output thresholding strategy to enable bispectrum phase es- 
timation at low light levels. 
In summary, we find that: 

(i) L3CCDs can be used for visibility estimation at any 
light level, though a small bias correction should be applied 
if a multiple threshold processing strategy is used on the 
L3CCD output. 

(ii) L3CCD visibility amplitude prediction at high 
(greater than about ten photons per pixel) and very low 
(less than 0.1 photons per pixel) photon rates is unbiased 
for all thresholding strategies. 
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Figure 4. Visibility phase error plotted against the number of fringe patterns recorded and processed when using an L3CCD as a 
detector. The solid black curves are an example of the phase error that may be obtained when using a perfect detector and closely 
matches Eq. [3] The dashed curves show the phase error with an L3CCD and an analogue thresholding strategy, the dotted curves for a 
PC thresholding strategy, the dot-dashed curve for a uniform thresholding strategy and the grey curves for a PP thresholding strategy. 
All plots are for a visibility amplitude of 0.1, and the mean light level is (a) 0.02 photons per pixel, (b) 0.2 photons per pixel, (c) 2 
photon per pixel and (d) 10 photons per pixel. We can see that the analogue thresholding strategy generally requires more photons to 
reduce the phase error to a given value. 



(iii) The use of multiple thresholding strategies will in- 
troduce a small bias for visibility amplitude at light levels 
between 0.1-20 photons per pixel. 

(iv) Since the mean gain is typically known to one per- 
cent accuracy, the error in visibility amplitude and phase 
estimation due to this uncertainty will be minimal. 

(v) At light levels greater than about five photons per 
pixel, the high-light-level bispectrum phase estimate can be 
obtained when using L3CCDs. At lower light levels down to 
about 0.1 photons per pixel, a bias term in the bispectrum 
phase becomes non-negligible, though this can be corrected 
by using the data to calculate the bias terms. Below this, a 
bias correction is still available when the raw L3CCD output 
is used, giving an unbiased estimate of the bispectrum phase. 

(vi) Using the raw L3CCD output with our adapted bis- 
pectrum bias correction, or a uniform thresholding strategy 
with the standard bias correction allows best bispectrum 
phase estimation. 

Our recommendation is that visibility amplitude and phase 
estimation at low light levels (less than 0.1 photons per pixel 
per readout) should be carried out using a single thresh- 



old processing strategy on the L3CCD output. At higher 
light levels, an analogue processing strategy should be used 
if there is ample signal-to-noise. If the SNR is low, using a 
multiple thresholding strategy with bias correction will help 
to improve the visibility estimate. Bispectrum phase esti- 
mation should be carried out using either the unthresholded 
L3CCD output or a uniform multiple thresholding strategy 
with the corresponding bias correction. 

Since L3CCDs provide the most accurate input estima- 
tion at low light levels using a single threshold, we recom- 
mend that, if possible, they are always used in this regime, 
increasing the frame rate if necessary to keep the number of 
photons per pixel low (< 0.1). We are currently developing 
a new controller for L3CCDs, allowing pixel rates of up to 
30 MHz, which will allow signal levels to be kept low in most 
interferometric applications. 
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Figure 5. Estimated bispectrum phase when using an L3CCD comparing the performance of different threshold processing strategies 
(represented by different curves) and estimation techniques (separate plots), (a) We sec that using the standard Wirnitzer correction, we 
are unable to estimate the bispectrum phase at low light levels when using analogue or PP thresholding strategies. However, by using 
our bias correction, we are able to use the analogue strategy to estimate the phase. We also see that the PC and Uniform thresholding 
strategies tend to overe stimate phase slightly at low light levels, (b) Using the bispectrum phase estimation technique developed by 
IPehlemann et alJ <1992l) we are unable to estimate phase accurately at light levels less than about 0.1 photons per pixel. The Poisson 
curve represents estimation when using a perfect detector, and the true bispectrum phase is 0.85 radians. 
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APPENDIX A: L3CCD UNIFORM 
THRESHOLDING STRATEGY 

iBasden et alJ ((2003) provide information about thresholding 
strategies which can be used to increase the SNR in a signal 
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from an L3CCD. For the purposes of this paper, we here in- 
troduce a new thresholding strategy which uses a knowledge 
of the mean light level to further improve the L3CCD signal 
SNR. The thresholds are placed with uniform separations, 
the separation being dependent on light level, chosen so that 
the mean photon flux is estimated correctly. 



Al Determining uniform threshold separations 

When using L3CCDs for interferometric signal detection, 
we are able to obtain additional information which will al- 
low us to develop a more accurate thresholding strategy. 
Since we make measurements of the interference signal over 
a fringe, we are able to estimate the mean light level within 
this fringe. We can then use this knowledge along with our 
knowledge of the mean gain to set the thresholds. Such a 
thresholding strategy is therefore dependent on the light 
level and so is only useful when an estimate of this can 
be obtained reliably in advance. We use a uniform threshold 
step size (with threshold boundaries being uniformly spaced) 
sized such that the mean thresholded output signal will be 
equal to the mean input (though we will not always thresh- 
old correctly for each event). We then interpret an L3CCD 
output signal falling between the n th and n + 1 th threshold 
as representing n photons. 

The probability of obtaining an output x electrons from 
the multiplication register of an L3CC D with a mean (Po is- 
son distributed) input, [i, is given by IBasden et aU l2003l: 



P(x) 



n=l 



cxp (— /j,)[i n x n 1 exp(—x/g) 
n!(n - l)\g n 



(Al) 



where g is the mean gain and /i is the mean light level. 

Thresholding this signal can improve the noise statis- 
tics. We can represent this by 



P(a) 



E 



P(x) 



(A2) 



where P(a) is the probability that a signal is put into the 
a th threshold, which has boundaries determined by f a . 



The expected output is then 



a = S ' aP(a). 



(A3) 



If we require a thresholding scheme which is able to pre- 
dict photon flux accurately without any additional correc- 
tions, we require the threshold boundaries to be dependent 
on light level. We equate the above equation to the mean 
light level, fi, using the threshold boundaries as our degree 
of freedom. If we wish to have uniformly spaced boundaries 
(f a = aT for threshold step size T), we can express this as 

°° r n 
exp /j, \ ^ J jj, 

M = ~g~ 2^\n\(n- 1)! 



E 



x=aT 

• E 

x=(a-l)T 



exp(-x/g) 



(A4) 



which should be satisfied by solving T for a given \i. 

A good approximation for the theoretical threshold step 



size at light levels greater than about 0.1 photons per pixel, 
T, is given by 



T ; 



9 



1 



2/i 3 / 4 



= 9 1 + 



2xV 4 



(A5) 



for gain g, and mean light level fi = x/g measured in ex- 
pected photons per pixel per readout estimated from an 
L3CCD output of x electrons. At light levels less than about 
0.1 photons per pixel it results in a larger threshold size than 
the ideal, though this is not problematic since at these light 
levels we are effectively using a large single threshold size, 
and input estimations are almost identical to those using a 
PC thresholding strategy. 

Using t his thresholding strat egy, the excess noise factor 
as defined by Basden et al is unity at low light levels, 

so the SNR scales as ^/n. This is because the step sizes here 
are large, and so we make correct predictions most of the 
time. As the light level increases, the e xcess noise factor 
tends towards \/2, the SNR scaling as y/n/2 (as with the 
analogue case), halving the effective quantum efficiency of 
the L3CCD. 

For cases where the mean light level can be estimated 
in advance of thresholding, this processing strategy leads 
to an improvement over the PC threshold strategy, being 
applicable at any light level, and results in an increased SNR 
over the analogue and PP strategies at low light levels, where 
the improvement is most welcome. It also has the advantage 
that no additional correction should be made after the event 
to preserve flux, since it is defined to estimate flux correctly. 



A2 Errors in threshold step size 

If the light level or gain are not known precisely, errors may 
be made when choosing the threshold step size. If the thresh- 
old step size is too small, we overestimate the photon count, 
while if it is too large, we underestimate. 

It is straightforward to show that the error in threshold 
step size depends on the error in gain and light level as: 



9 
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(A6) 



where A represents the uncertainty in a measurement, and 
the second equality is valid if the mean light level fi is calcu- 
lated from m averaged raw L3CCD outputs, each containing 
a mean of n photons with a SNR of yn/2 due to the Poisson 
and multiplication processes. We can u sually expect to know 
the mean gain to about one percent iMackav et aljEoOlT) . 
and so this error is the limiting factor when m 3 n> 10000, 
for example n = 0.0005 when m = 256. We can therefore 
usually expect to be able to choose our threshold step size 
with an accuracy equal to that in the mean gain. 

At high light levels, the fractional error in prediction is 
proportional to the fractional error in threshold step size. 
At low light levels, there is less of a dependence since the 
thresholds are widely spaced, and signals usually fall within 
the first threshold. 
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APPENDIX B: MULTIPLE THRESHOLD 
VISIBILITY AMPLITUDE BIAS 

If we use multiple thresholding on the L3CCD output we 
will introduce a systematic bias into our visibility ampli- 
tude estimations, though we can counter for this effect since 
we can determine the bias. The reason for this bias is now 
discussed. 

Bl PP thresholding 

When we use a PP thresholding strategy, to estimate the 
flux correctly we m ust apply a correction according to 
feasden et alJl2003ft 

h/ It « [1 + 0.7 exp(- Jc/3)] . (Bl) 

with subscripts e and r corresponding to estimated and real 
values respectively. Since this correction depends on the light 
level, we are therefore biasing our visibility amplitude esti- 
mation since fringe minima will have a different weighting 
(requiring a different correction) than fringe maxima. Our 
estimate of visibility amplitude will therefore be biased. 

We are able to verify this by considering the definition of 
visibility amplitude, V = 7 + +J -~ where I± is the intensity at 
the maximum and mini mum of the fringe. Since 1 < J c // r < 
1.7 feasden et al .120031) . we can further approximate Eq. lBll 
to 

/ e «ir [1 + 0.7 exp(-J r /3)]. (B2) 

Therefore, using our estimated flux with I r ± = /i(l ± V T ) 
for mean light level fi, and inserting into the definition of 
visibility amplitude, we will obtain an estimate (V e ) of the 
true visibility V T ) equal to: 

v „ 2K + 0.7 [U exp (-/!/+ /3) - J- exp (-/i/-/3)] , . 
c ~ 2 + 0.7[/+exp(-M/+/3)+/_exp(- M /_/3)] 1 ' 

where f± = 1 ± V r . We find that this agrees very closely with 
our calculation using the probability distributions, Eq. [5] 

B2 Uniform thresholding 

If we use a uniform thresholding strategy, then since we are 
able to measure the mean intensity (fi photons per pixel) 
of an interference fringe we use this information to set the 
threshold step size. However, for visibility V, the fringe in- 
tensity will actually vary from a maximum of pi(l + V) to 
a minimum of — V). At these points our threshold step 
sizes will be wrong. At a minimum of intensity, the threshold 
step size chosen from the mean light level will be too small 
leading to an overestimation of signal, while at intensity 
maxima the threshold step size will be too large resulting 
in an underestimation. The combined effect is to underesti- 
mate the visibility. 

B3 Low and high light levels 

Visibility amplitude bias correction is not be required for 
mean light levels below about 0.1 photons per pixel, and 
above about 10 photons per pixel. At the lowest light levels, 
we will generally only be using the first one or two thresh- 
olds, and every photon event will have the same chance of 



falling into these. At higher light levels, we find that the 
bias in our visibility amplitude estimation will tend towards 
zero, leading to correctly estimated visibility amplitude. 



APPENDIX C: BISPECTRUM PHASE BIAS 
ESTIMATION 

CI Unbiased bispectrum phase estimation: 
Method I 

IPehlemann et al] lll992l) have studied photon bias compen- 
sation of both power spectrum and bispectrum estimation, 
and derive a method for correcting this. They assume that 
an image is made from individual photons, each of which 
may have a different peak intensity, ctj. We can relate this 
to the L3CCD output, using the L3CCD output probability 
distribution. The total expectation value E[D^ 2 '(u)] of the 
power spectrum of a data set is given by 

E[D {2) (u)] = cps,o < J {2) (it) > +c PS ,i (CI) 

where cps.o = N 2 ~a 2 and cps,i = Na 2 are scalar coefficients, 
and < J^ 2 \u) > is the average power spectrum of the nor- 
malized high-light-level (unbiased) fringe. We find that this 
expression can be applied to the L3CCD raw output since 
it has this form. However, when using other L3CCD thresh- 
olding strategies, this expression is not strictly true, though 
we do indeed get a flat noise background, cps.i. 

The t otal expectation fo r the b ispectrum, E[D^ (it, v)] 
is given bv IPehlemann et al] il992T) as 

E[D (3) (u, v)] = cbs,o < J (3) (it, v) > +C BS ,2 

+ C BS ,i(< J {2) (u)> + < J {2) (V)> 

+ < J (2) (-u-i;) >) (C2) 

where cps.o = N^a 3 , cbs,i = N 2 aa 2 , cbs,2 = Na 3 are scalar 
coefficients and < (u, v) > is the desired average bispec- 
trum of the normalized high-light-level (unbiased) fringe. 

Since we know that the probability distri bution for the 
L3CC D output given a single photon input feasden et al] 
2003) is P(x) = cT 1 exp{-x/g), we can calculate the coeffi- 
cients as: 

1U 2 2 
cps,o — JV g 

c PS ,i = 2Ng 2 

"at 3 3 

cbs.o = N g 
cbs, i = 2N 2 g 3 
cbs,2 = 6Ng 3 

where g is the mean gain. When using the analogue thresh- 
olding strategy, we simply divide the L3CCD output by g. 
This allows us to write 

N S < J (3) (it, v) > = E[D {3) (u, v)] + 67V 

- 2 {E[D (2) (it)] + E[D (2) («)] 

+E[D (2) (-u~v)]) (C3) 

from which we can obtain an unbiased estimate of the bis- 
pectrum phase. 
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C2 Unbiased bispectrum phase estimation: 
Method II 

If we are not easily able to calculate the probability distri- 
bution for a, we can estimate the required correction terms 
using our data. To do this, we evaluate the bispectrum at 
many frequencies chosen so that we can obtain the required 
coefficients one at a time. 

(i) By choosing two frequencies, /i and f% such that 
fi 1/2 ^ u,v,u + v where it, v and u + v are the frequencies 
at which the peaks appear in the power spectrum, we can 
see that E[D^ (/1, /2)] = cbs,2, and so obtain an estimate 
for c B s,2- 

(ii) Provided u 7^ v and u 7^ 2v we evaluate the bis- 
pectrum at /1 = /2 = u,v,u + v. At these frequencies, 
< J^ 3 ' (/1, ^2) >= and so the bispectrum here evaluates to 

CBS.2 + 2C SS ,1 < J (2) (/l) > +CBSA < J (2) (2/l) >• 

(iii) Using our power spectrum (Eq. lCH we can estimate 
cps,o < j( 2 '(u) >= E[D^ (u)] — cps,i, obtaining cps,i from 
the flat noise background. 

(iv) We are now able to estimate cbs,i/cps,o from the 
previous two steps: 

CB5.1 _ E[D (Z) (u,u)] - C B S,2 , C4 > 

cps,o ~ 2E[DW(u)] -3cps,i 

(v) The unbiased bispectrum phase can then be obtained 
at fi — u and f% = v, since this is independent of css.o- 

This paper has been typeset from a TfjSC/ I^TjjX file prepared 
by the author. 



